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Abstract. The aim of this paper is to propose methods that enable us to 
build new numerical schemes, which preserve the Lie symmetries of the original 
differential equations. To this purpose, the compound Burgers-Korteweg-de 
Vries (CBKDV) equation is considered. The particular case of the Burgers 
equation is taken as a numerical example, and the resulting semi-invariant 
scheme is exposed. 



1. Introduction. Finite difference equations used to approximate the solutions of 
a differential equation generally do not respect the symmetries of the original equa- 
tion, and can lead to inaccurate numerical results. Usually, specific equations are 
considered, for which the authors build a scheme preserving the symmetries of the 
original differential equation. Yet, it is more interesting to directly consider a class 
of differential equations, in order to obtain more general results. 

Using the work of Yanenko [4] and Shokin [3] , who applied the Lie group theory to 
finite difference equations by means of the differential approximation, in conjunction 
with the approach of Ames et al. [8], we generalize results developed in pQ, [2], and 
expose the invariance condition for a differential approximation of the CBKDV 
equation. A numerical example, in the particular case of the Burgers equation, is 
developed. 

2. Notion of symmetry. 
2.1. Definitions. 

Definition 2.1. A r-parameter Lie group is a r-dimensional smooth manifold GV, 
which has the group properties, such that the group operation of multiplication and 
inversion are smooth maps. 
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Especially, a Lie group G r is defined as a group of continuous transformations which 
act on an open subset of the Euclidean space R fc of variables, which change under 
the action of G r . We presently concentrate on a local group, the transformations 
of which are close to the identity transformation. 

Definition 2.2. A r-parameter Lie group G r is a group of point transformations, 
which acts on X x U, the space of the independent variables and the dependent 
ones: 

G r = {x* = <f>i(x,u, a); u* = ipj(x,u,a), i = l,...,m; j = 1, . . . , n; a = (oi, . . . , a r )} (1) 
where x £ X C K m and u £ U C R n . 

G r locally satisfies the group axioms: existence of an identity element, associativity, 
inversibility, closure under the binary composition operation. The transformation 
corresponding to a zero parameter is the identity transformation. 
Expand the transformations by means of a Taylor series at the zero value of the 
parameter a: 

Xf =Xi + a a - — +0(a a ), a = l,...,r 
aa a \a=o 



u* =Uj+a a ^p-\ +0{a 2 a ), a= l,...,r (2) 
aa a \a—o 



The derivatives of <\>i and ipj with respect to the parameter a a are smooth functions, 
called infinitesimals of the group G r . Denote by £f and rf? the infinitesimals of G r . 
The point transformation group G r can be represented by means of the operator 
L Q : 

d_ . „, , d 
' dxi 



£f(x,u)— \-r)"(x,u) , i = l,...,m; j = l,...,n; a = l,...,r (3) 



The operators L Q , a = 1, . . . , r are called the infinitesimal operators of G r . 
{L Q , a = 1, . . . ,r} represents the set of tangent vectors to the manifold G r , when 
the zero value is assigned to the parameter a. The set is a basis of the Lie-algebra 
of the infinitesimal operators of G r , the dimension of which is the same as the one 
of the Lie group G r . 

The knowledge of the L a enables us to determine the point transformations of the 
group G r by solving the equations: 

dx* du* 

Tr L= ^( x ^ u *)' = n?(x*,u*), i = 1, .. . ,m; j = 1, .. . ,n; a = 1, . .. ,r (4) 

aa a oa a 

in conjunction with the initial conditions: 

(5) 



Example 2.3. The Galilean transformation group 

The Galilean transformations correspond to time dependent translations of a refer- 
ence frame : 

G = {T : (x, t, u) i — ► (x*, t*, it*) = (x + et, t,u + e)} 
where e is the translation constant velocity. 

T and its inverse function are continuous. The infinitesimals functions of the group 
axe = 0,1) 
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2.2. Symmetry properties of differential equations. The notion of symmetry 
is a tool for generating new solutions of differential equations. Let us review the 
main aspects of the application of the Lie group theory to differential equations. 
Consider a system of Z'^-order differential equations: 

F x {x,uM k:L \^ klM K---,ii ( - kl '" k ^) = 0, X = l,...,g (6) 

Denote by u ( fe i - fc p) the vector, the components of which are partial derivatives of 
order p, namely, nf 1 '"^ = g^^d^T = • • • > n and k\, . . . ,k p £ {1, . . . , m}. 
Denote by x = {x\, . . . , x m ) the independent variables, u = [u\, . . . , u n ) the depen- 
dent variables, and (xk ± ■ ■ ■ Xk p ) a set of elements of the independent variables. 
Equation J6|) is a subset of X x U' 1 ' , a prolongation of the space X xU to the space 
of the partial derivatives of u with respect to x up to order I. X x [/W, which is a 
smooth manifold, is called the l-th order jet space of X x U. In order to take into 
account the derivative terms involved in the differential equation, the action of G r 
on X x U needs to be prolonged to the space of the derivatives of the dependent 
variables. 

Denote by Gr a r-parameter Lie group of point transformations acting on an open 
subset M^ 1 ' of the Z-th order jet space X x J/W of the independent variables x, 
dependent variables u and the partial derivatives of u with respect to x. 
The ^-prolongation operator of G r is: 

i = 1, . . . , m; j = 1, . . . , n; a = 1, , r. 

The infinitesimal functions £f , a"''^ and a -J>( kl — ko ' are gi ven by: 

w _ a< fc I v a _ d( Pi I o-?»( fc i) - V' ( 8 ) 



p g , 0=2 , 

' Dx k iti 9xidx kl . . . dx ko _ 1 T>x ko ' 



, V d J\ du.; d 
where: = h > — 

fit ozfc ^ a^fc ouj 

Theorem 2.4. TTie system of I th -order differential equations is invariant under the 
group G r if and only if it satisfies the following infinitesimal invariance criterion: 

Z ( ^T X \ =0, a = 1,. . . ,r; A = 1,. . . ,q (9) 



3. Lie group of a differential approximation. The symmetry group analysis 
of the differential approximation uses the techniques of the Lie group theory applied 
to differential equations. The differential approximations involve step size variable 
which change under the group action. 

The technique of symmetry analysis is not directly applied to the finite difference 
schemes. It is based on the differential approximation, which describes approxi- 
mately the numerical solution behaviour at a reference point of the mesh. Thus the 
concept of differential approximation is a local object, which can not systematically 
detect a mesh change. Despite the local behaviour of the symmetry group analysis, 
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the invariant method based on the differential approximation has enabled one to 
establish interesting properties of symmetry group of finite difference schemes. 
The finite difference scheme, which approximates the differential system |(6|), can be 
written as: 

A x (x,u, h,Tu) = 0, X = l,...,q (10) 

where h — (h 1 ,h 2 , ■ ■ ■ , h m ) denotes the space step vector, and T = (T 1; T 2 , . . . , T m ) 
the shift-operator along the axis of the independent variables, defined by: 

Ti[u](xi,X2, ■ ■ . , Xi-i,Xi, Xi+l, ■ ■■ ,X m ) = u(xi,X2, ■ ■ . ,Xi-l,Xi + hi,X i+1 , ... ,x m ). (11) 

Definition 3.1. The differential system: 

V X {x,u,u^\...,u^- k ^) = T X {x,u,u^\...,u^- k ^) 

s m 

+E X»^0w (fci) ,- ■ .,» (fcl - k "^ ) ), 

0=1 i=l 

X = l,...,q; I' = mctX(x,i)l' A,i (12) 

is called the s th -ovdev differential approximation of the finite difference scheme 
iflOj) . In the specific case s = 1, the above system is called the first differential 
approximation. 

The differential system ((3J) is obtained from the algebraic system lfT0|) by applying 
Taylor series expansion to the components of Tu about the point x — {x\, . . . , x m ) 
and truncating the expansion to a given finite order. 

Denote by G' r a group of transformations acting on an open subset M' of X x U x H 
the space of the independent variables, the dependent variables and the step size 
variables : 

G' r = {x* = <j>i(x,u, a); u* = ipj(x,u,a);h* = ipi(x,u,h,a), i = 1, . . . , m; j = 1, . . . , n} (13) 
by Jj a the basis infinitesimal operator of G' r : 

= ^(x,u)^+ V f(x,u)^-+^(x,u,h)^-, a = l,...,r (14) 
axi J auj ahi 

where Cf = ttA , a=l,...,r 

and by Gr ^ a group of transformation acting on an open subset ' of the space of 
the independent variables, the dependent variables and the step size variables and 

the partial derivatives involved in the differential system . The ^''^-prolongation 

operator of G' r , Lq ' can be written as: 

Theorem 3.2. The differential approximation (E) is invariant under the group 
Gr ' if and only if 

£2' ) -p A ((a;,u,u( fsi ),...,u( fcl -V))| ^ = 0, a = 1, . . . ,r; A = 1, . . . , q (16) 

s m 

or [lSV+lJI'^^C 1 ^^)]! a= = °> a=l,...,r; X = l,...,q (17) 

/3 = 1 i=l 

4. Determination of the infinitesimal functions. The calculation of Lie groups 
of differential equations with pencil and paper is tedious and may induce errors. The 
size of related equations increases with the number of the symmetry variables, and 
the order of the differential equations. A large amount of packages have been cre- 
ated using software programs with symbolic manipulations, such as Mathematica, 
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MACSYMA, Maple, REDUCE, AXIOM, MuPAD. Schwarz [19] wrote algorithms for 
REDUCE and AXIOM computer algebra systems, Vu and Carminati [16j worked on 
DESOLVE, a Maple program, Herod [17] and Baumann [18] developed Mathematica 

programs. 

The authors have implemented a Mathematica package [1], for the determination of 
the Lie group of the differential approximation of one dimensional model equations. 
Theorems 12.41 and I3.2[ respectively give the algorithmic procedure for the determi- 
nation of Lie group of any differential system and differential approximation. 
The theorem infinitesimal invariance criteria involved the independent variables x, 
the dependent ones u, products of the partial derivatives of u with respect to x, the 
unknown infinitesimal functions £f , and 77", i — 1, . . . , m; j = 1, . . . , n and their 
partial derivatives with respect to x and u. 

The partial derivatives of the infinitesimal functions with respect to x and u come 
from the coefficients of prolonged Lie algebra vector field (JHJ • 

Equations ([9]) and lfT7|) are simplified by means of the conditions ^ and ((3]) . This 
simplification manipulation eliminates some derivatives of u. These equations are 
then solved algebrically with respect to the partial derivatives of the dependent 
variables, handled as independent variables. Denote by w the vector, the compo- 
nents of which are these variables. Since the whole equation holds for all the w 
components, each coefficient in front of the products of the w components has to be 
zero. This leads to a linear overdetermined system of partial differential equation, 
with respect to the infinitesimal functions, called the determining equations of the 
Lie group of the related differential system . 

The solution of the overdetermined system can be found either by using elementary 
methods of the theory of linear partial differential equations or by using a poly- 
nomial form for the infinitesimals. The last technique provides a linear system of 
algebraic equations with respect to the polynomial coefficients. But it can not find 
infinitesimals with transcendal functions. 

The resolution of the determining equations yields explicitly the expression of £f , 
rjf, a = 1, . . . , r, i = 1, . . . , to, j = 1, . . . , n. Then relations and (J5j) provide the 
calculation of group transformations from the infinitesimal expression. 

5. S ymmetries of the CBKDV equation. The invariance condition of theorem 
I2.4[ which enables us to obtain the expression of the infinitesimal operators, can be 
written as: 

LL 2 V| =0, a = l,...,r; (18) 

Equation l[T8|) provides an overdetermined system of linear partial differential equa- 
tions for the infinitesimal functions. 

Solving these equations, by means of a symbolic calculus tool (Mathematica), pro- 
vides the expression of the infinitesimal functions of a two-parameter group: 

& = ^0 , & = b , rf = 

where ao, bo are constants. 

The two-dimensional Lie algebra of the group G is generated by the following op- 
erators: 

T 9 T 9 

dx dt 

which respectively correspond to: 
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• the space translation : (x,t,u, u) i — > (x + xo,t, u, v)\ 

• the time translation : (x,t,u, u) \ — > (x,t + t ,u,u); 

xq, to are constants. 

It is important to note that, in the case of particular cases of the CBKDV equation, 
when, for instance, some of the coefficients are equal to zero, more symmetries can 
be obtained. 

6. The specific case of the Burgers equation. Yanenko [4] and Shokin [3] 
symmetry analysis is applied to finite difference schemes for solving the Burgers 
equation. Thus, the symmetries of the Burgers equation, which are finally broken 
by the finite difference discretization, are determined. The techniques exposed in 
[3] and [4] enable one to construct differential approximations, which preserve the 
symmetries of the original differential system. We call the related finite difference 
scheme a semi-invariant scheme, in so far as the invariance condition is weaker 
than the one of the other invariance methods, defined as direct invariance methods. 
Indeed, the approach in [3 J and |4J does not deal with the invariance of the algebraic 
equations, which govern the mesh evolution. 

The differential equation has provided important characteristics for numerical schemes, 
in the study of numerical stability, dissipation and dispersive property. In [3j, (4] 
and |8j, the differential approximation has been revealed as a practical and recevable 
tool for symmetry analysis of finite difference scheme. 

A comparison is made between the numerical solutions of the Burgers equation for 
some standard schemes and the semi-invariant one. 

6.1. Symmetries of the Burgers equation. The Burgers equation, which is a 
particular case of the CBKDV equation, with a = 1, j3 = 0, s = 0, fi = —is, can be 
written as: 

T(x, t, u, v, u x ,ut, Mm) = tit + u u x — v tin = (19) 
where v > is the dynamic viscosity. 

Let us denote by G a group of continuous transformations of the Burgers equation 
acting on an open subset M of the space the independent variables (x,t), the depen- 
dent variable u, and the viscosity v. The viscosity is taken as a symmetry variable 
in order to enable us to take into account variations of the Reynolds number. 
The six-dimensional Lie algebra of the group G is generated by the following oper- 
ators: 

d d a a a 

Li = — , Li2 = — , L3 = x h 2t u — 

dx dt dx dt du 

a 9 a . a a a a a a , , 

L 4 = xt ht h(-ut + x) — , L 5 =t 1 , L 6 = -t hii Yv — (20) 

dx dt du dx du dt du du 

which respectively correspond to: 

• the space translation : (x, t, u, v) \ — > (x + ei, t, u, v)\ 

• the time translation : (x, t, u, u) 1 — ► (x, t + e 2 ,u, u); 

• the dilatation : (x,t,u, v) 1 — ► (e 3 x, e§t, 1 u, u); 

• the projective transformation : {x,t,u,v)\ — > ( 1 _^ t , 1 _' e4t , 3^4 + "(1 - ^t), u^j ; 

• the Galilean transformation : (x,t,u, v) 1 — > (x + <e 5 t,t, u + es,u); 

• the dilatation : (x,t,u, v) \ — > (x, <=g 1 t, e 6 u, e 6 u). 

(ei)i=i,...,6 are constants. 

As expected, the Burgers equation, as a particular form of the CBKDV equation, 
admits more symmetries, yieldind a richer Lie algebra. 
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6.2. Symmetries of first differential approximations. Denote by h the mesh 
size, t the time step, N x the number of mesh points, N t the number of time steps, 
and it™, i 6 {0, . . . , Nt}, n S {0, . . . , N x } the discrete approximation of u(ih, tit). 
In order to shorten the size of the finite difference scheme expressions, we use the 
following notations introduced by Hildebrand in |20| : 

U n 1 - 1 U" j + U n _ i 

s«)= - MO = 2 2 1 5 

U n , 1 — U n 1Z n — U" ! 

5 K ) = t . 5 K ) = r 1 s "i = 

The Burgers equation can be discretized by means of: 

• the FTCS (forward-time and centered-space) scheme: 

T h V 2 '« h 2 1 



• the Lax-Wendroff scheme: 

71+1 7) f 



h y 2 n h 2 



where: 



Af= _ - [ E l <5+( ^ r _ £ ,-| u n 5 - ( ^ ) «i_^ir^„ 

• 2 /i 2 L 1 y 2 l% y 2 ll J 2 U 4 



2/i 3 L ! ' 1 ' 1 v J 2 L h 3 v 2 

the Crank-Nicolson scheme: 



1 ii 

The linear stability properties and the related orders of approximation are: 

• the FTCS scheme: S<\, CFL < 1, CFL 2 < 2S; 0(r,h 2 ) 

• the Lax-Wendroff scheme: S* < §, CFL < 1; 0{r 2 ,h 2 ) 

• the Crank-Nicolson scheme: unconditional stability; 0(r 2 ,h 2 ) 
where CFL = Sfc, 5 = $ and 5* = {v + ^^r^) pr- 

Consider function of the time step r, and of the mesh size h, expand it at 

a given order by means of its Taylor series, and neglect the o(r Q ) and o(h^) terms, 
where a and f3 depend on the order of the schemes. This yields the differential 
representation of the finite difference equation. 
The following differential representations are obtained: 

• for the FTCS scheme: 

1 o r h 2 2 uh 2 

Ut + -(« )x - V Uxx + — <?2 + T^( u )xxx —Uxxxx = 

^ _ _L _L 

• for the Lax-Wendroff scheme: 

1, 2 , -r 2 h 2 , 2 s vh 2 
Ut + -(it ja; — f Usa + —S3 + T^( u )»« —Uxxxx = 

• for the Crank-Nicolson scheme: 

1 2 n 2/93 1/2 \ ^ / \ \ ^ 2 /2\ 
Ut + -(« )a; — i/Mm + r (— + -(g 1 + ug 2 )x - -(92)xx) + — [u ) X xx — %« = 

where 9l = -(^) x + uu X x, g 2 = (- giu) x + u( gi ) xx , g 3 = ( - g 2 u - g 2 ) x + v{g 2 ) xx 

Denote by G' the group of transformations of a first differential approximation 
acting on an open subset M' of the space of the independent variables (x,t) and 
the dependent variable u, the step size variables (h, r) and the viscosity v. 
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The /'^-prolongation of G' can be written as: 



i' 



9x (9t 9m ^- ' a, (fci...fc p ) g/j g r ,9^ 



p=l 



where /' h as be en defined in definition 13.11 

Theorem 13.21 enables us to obtain the necessary and sufficient condition of invari- 
ance of the first differential approximation V: 



% l ) v\ = (22) 



Theorem 13.21 is applied to the differential representations of the above schemes. 
The resolution of the determining equations of each first differential approximation 
yields the 4-parameter group (see [1]): 

£%=a + bx, & = c+(2b-d)t, yf = (-6 + d) u (23) 
£f = b ft, C2 = (26 - d) r, 6» a = ei/ 



The 4-dimensional Lie algebra of G' is generated by: 

9 T 9 T , 9 9 8,8 8 
Li = — , L 2 = — , I/3 = x h 2t u \-h h2r — 

dx dt dx dt du dh dr 

1 d d d d 

I/4 = — t \-u t Yv — (24) 

dt du dr du 



These operators are respectively related to: 

• the space translation : (x,t,u, h,r, v) 1 — > (x + ei,t,u, h,r, u); 

• the time translation : (x,t,u, h,r, v) — > (x,t + e2,u, ft, r, v); 

• the dilatation : (x,t,u,h,T,u) 1 — > (e^x, e§i, €3 1 u, ezh, e\r, v); 

• the dilatation: (x,t,u,h,r,u) 1 — > (z, 1 t, <E4M, ft, 1 t, <E4^); 

where (61)1=1^.^4 are constants. 

The above finite difference equations are preserved by the space translation, the 
time translation and both dilatations. 

Approximating the Burgers equation by the above finite difference equations results 
in the loss of the projective and Galilean transformations. 

6.3. A semi-invariant scheme. Two analogous methods provide a direct sym- 
metry analysis of finite difference schemes and can lead to the definition of adapted 
evolutionary meshes, whose geometrical structure is preserved by the entire group. 
The first direct method has been introduced by Dorodnitsyn [IT] and is based on 
Lie algebra techniques, using the infinitesimal operators. The second method has 
been introduced by Olver [5] and is based on the theory of the Cartan moving frame. 
The method proposed by Yanenko [4] and Shokin [3] consists of a symmetry study 
of the differential approximation. Although the last method is not fully exact, the 
numerical results in [3] and [8] has proved its effectiveness. 
The scheme proposed below is associated to an uniform orthogonal mesh. 
We propose to approximate the Burgers equation by the finite difference scheme: 

< +1 ~ «" + 1 ( _ Ajg.) /«V _ 1 , g 2 _ £\ _ h ( Q n s + _ Q n 5 -\ K=Q (25) 

r ft v 6 JK 2 Jt ft 2 v 12 J 1 V »+± »-i J 1 y ' 

where O™ = Vt(xi,t n ,uf) is defined next so that the related differential represen- 
tation is preserved by the symmetries of the Burgers equation. The scheme has 



TOWARDS NEW SCHEMES: A LIE-GROUP APPROACHOF THE CBKDV AND ITS DERIVED EQUATIONS 



second-order accuracy in space and first-order accuracy in time. The derivatives 
(u 2 ) x and u xx are approximated by fourth order accuracy difference expressions: 

-^)(«?)=(«--^«B-)r+o(h B ), (S-i2^) ( < )= ^-S Mfa) " + ° (?i6) (26) 



The truncation error of the difference scheme l|25]l can be written as: 

e = "-uu - h 2 (nu x ) + 0(r 2 ) + 0{h A ) 
2 V / x 

u u is replaced by an expression involving partial derivatives with respect to x, by 
using the Burgers equation. Replacing the obtained expression in the truncation 
error leads to: 

/ \ VT ( \ VT ,U 2 . V 2 T , o , A . 

e = (Cu x ^-—( y uu xx ) x ~—{—) xxx + —u xxxx + 0(r 2 ) + 0(h 4 ) 

where C = ^u 2 — h 2 Q. 

It is convenient for the calculation of C that the truncation error is reduced to: 

e = (Cu x ^+0(r 2 ) + 0(h 4 ) 

The related finite difference scheme is the following first order accuracy in time and 
second order accuracy in space: 

7/ n + 1 — li n 1 //A 3 7/ 2 1 A 4 / \ 

h !k + -(ns- ) n - u^-(s 2 - -)«) - h(n\ ,f+ - n n 

t h y 6 JK 2 >% h 2K 12' y 1 ' V '+i ) 1 

2 V h 2 K h 2 y '~rJ 2 h A 1 2 K 2 U K ' 

and the differential approximation can be written as: 

V(x,t,u,v,u x ,ut,u xx ) = u t + u u x — v Uxx + (Cu x ) x = (28) 

The von Neumann stability analysis of scheme (|27| under a linearized form provides 
the following necessary conditions for S, CFL and fl T = fir: 

CFL 2 — 2S — 2Q T < 0, < (4S)/3 - 2S 2 + Q T < 1/2 (29) 

If £1 is sufficiently close to zero, these conditions become then sufficient for the 
linear formulation. 

6.4. Numerical application. The artificial viscosity term (Cu x ) x is build so as 
to preserve the symmetries of the Burgers equation. C is function of (x, t, u, r, h). 

The Burgers equation is solved numerically for the semi-invariant scheme, and the 
classical schemes mentioned above. The numerical solutions are determined in the 
frame (Fl), with an uniform and orthogonal mesh, in the frame (F2) and (F3), 
which are respectively obtained with the Galilean transformation [x, t, u, v) i — > 
(x + 0.25 t, t, u + 0.25, z/) and a higher translation velocity (x,t,u,v) i — ► (x + 
0.5 t, t, u + 0.5, v). In (F2) and (F3) the orthogonality of the mesh is not preserved. 
The parameters, which control the numerical stability, i.e. the CFL number and the 
cell Reynolds number, are preserved by the Galilean transformation. Any change of 
the frame does not modify their value. These parameters satisfy the linear stability 
conditions, which are sufficient for all the numerical applications realized below. 

Figures Q] display the time evolution, of the L 2 norm of the error, for the above 
numerical schemes, in the different frames. A shock wave type solution is initially 
introduced. The CFL and the cell Reynolds number are respectively equal to 4 
and 0.5. Non-invariant schemes are sensitive to the change of frame. Thus, the 
numerical error of the FTCS scheme oscillates after the application of a translation 
with a velocity of 0.25. It increases sharply after a translation with a velocity of 0.5. 



10 



EMMA HOARAU, CLAIRE DAVID 




Figure 1. Evolution in time of the L norm of the error in (Fl), (F2), (F3) when 
Re h = 4.0, CFL = 0.5 

The linear numerical stability is satisfied, though the discretization parameters of 
the CFL and the cell Reynolds numbers are situated on the boundary between the 
region where the solutions rapidly blow up and the one where the solutions remain 
bounded for certain time steps. The introduction of an additional error due to the 
loss of the Galilean transformation leads to the numerical solution blow-up. 
We notice the dependence of numerical error, for the semi-invariant scheme, upon 
the frame. The change of the frame causes the apparition of slight oscillations. 
Figures [2] display the time evolution, of the L 2 norm of the error, for the above 
numerical schemes, in the same frames. This time, a wavy solution is initially 
introduced. The CFL number and the cell Reynolds number are respectively equal 
to 0.2 and 0.06. We clearly notice the dependence of numerical error, for the four 
numerical scheme, upon the frame. 

7. Conclusion. Starting from the general CBKDV equation, we have set methods 
that enable us to build new numerical schemes. The application to the Burgers 
equation can be seen as a restricted case of the CBKDV equation. 
The invariance method based on the differential approximation does not enable us 
to build schemes, which are rigorously invariant under the action of the group of the 
original equations. Like in the case of classical finite difference methods for aerody- 
namics, the method based on the differential approximation breaks the symmetries, 
which act on the independent variables, and transform the geometrical properties 
of the mesh. This symmetry action is observed for the Galilean transformation and 
has been applied for the numerical resolution of the Burgers equation. Under the 
action of the Galilean transformation, a non-invariant discretization leads to the 
loss of the mesh orthogonality. 
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FIGURE 2. Evolution in time of the L 2 norm of the error in (Fl), (F2), (F3) when 
Re h =, CFL = 

If, on the one hand, the method based on the differential approximation would 
provide the exact invariance if the higher order terms of the numerical error were 
considered, it would yield, on the other hand, a scheme too complicated to be 
implemented. In order to preserve rigorously the symmetries, we need to use direct 
methods of symmetry analysis of finite difference schemes, which take into account 
the transformation of all the finite difference variables and preserve the geometrical 
properties of the mesh. 
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